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Abstract 

Distributions exhibiting fat tails occur frequently in many differ- 
ent areas of science. A dynamical reason for fat tails can be a so- 
called superstatistics, where one has a superposition of local Gaussians 
whose variance fluctuates on a rather large spatio-temporal scale. Af- 
ter briefly reviewing this concept, we explore in more detail a class 
of superstatistics that hasn't been subject of many investigations so 
far, namely superstatistics for which a suitable power (3^ of the local 
inverse temperature /3 is x'^'distributed. We show that r/ > leads 
to power law distributions, while r/ < leads to stretched exponen- 
tials. The special case rj = 1 corresponds to Tsallis statistics and the 
special case r] = —1 to exponential statistics of the square root of 
energy. Possible applications for granular media and hydrodynamic 
turbulence are discussed. 



1 Introduction 



Nonextensive statistical mechanics ^1 |2l El E] was originally developed as an 
equilibrium formalism, but most physical applications of this formalism ac- 
tually occur for typical nonequilibrium situations. Sometimes these nonequi- 
librium situations are described by a fluctuating parameter /3, which may, 
for example, be the inverse temperature. Alternatively, jS may be an effec- 
tive friction constant, a changing mass parameter, a changing amplitude of 
Gaussian white noise, a fluctuating local energy dissipation or simply a local 
inverse variance parameter extracted from a time series. The fluctuations of (3 
induce a superposition of different statistics on different time scales, in short a 
superstatistics 013 13 li 011011111 1121 113 III 113121 113 III IIS 1201 IH^ 
The stationary probability distributions of superstatistical systems typically 
exhibit much broader tails than a Gaussian distribution. These tails can 
decay e.g. with a power law, or as a stretched exponential, or in an even 
more complicated way 8j. Which type of tails are produced depends on the 
probability distribution /(/?) of the parameter f3. Recent applications of the 
superstatistics concept include a variety of physical systems. Examples are 
Lagrangian 1211231231211 and Eulerian turbulence [231231301, defect tur- 
bulence atmospheric turbulence cosmic ray statistics solar 
flares [23 , solar wind statistics [201 , networks jSZlEHli random matrix theory 
jnOj, and mathematical finance [131111 1121 ■ 

If (3 is distributed according to a particular probability distribution, the 
X^-distribution, then the corresponding marginal stationary distributions of 
the superstatistical system obtained by integrating over all /? are given by 
the generalized canonical distributions of nonextensive statistical mechanics 
[II 13 13 HI- For other distributions of the intensive parameter /3, one ends up 
with more complicated statistics. 

In this paper, after briefly reviewing the superstatistics concept, we ex- 
plore a rather general case which may be of relevance to many practical 
applications. We consider the case of a superstatistics where i.e. P to 
some power t], is x^-distributed, where rj is some arbitrary parameter. The 
case ?7 = 1 is fully understood: It leads to Tsallis statistics and asymptotic 
power-law decay of the marginal distributions obtained by integrating over 
all p. However, the other values of rj are interesting as well, and will be ex- 
plored in more detail here. For general rj > Owe obtain asymptotic power law 
decay, though the resulting statistics is slightly different from Tsallis statis- 
tics (only 7] = 1 leads exactly to Tsallis statistics). For rj < one obtains 
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tails that asymptotically decay as stretched exponentials. The special case 
?7 = — 1 corresponds to exponential tails of the square root of energy. We 
will provide some arguments (based on the ordinary Central Limit Theorem) 
why nonequilibrium systems with many degrees of freedom often lead to one 
of the superstatistics described above. 

2 Various types of superstatistics 

It is well known that for the canonical ensemble the probability to observe a 
state with energy E is given by 

m = ^P(^)e-'^. (1) 

g-/3E jg ^Yic Boltzmann factor, p{E) is the density of states and Z{P) is 
the normalization constant of p{E)e~f^^ . For superstatistical systems, one 
generalizes this approach by assuming that /3 is a random variable as well. 
Indeed, a driven nonequilibrium system is often inhomogeneous and consist 
of many spatial cells with different values of j3 in each cell. The cell size 
is effectively given by the correlation length of the continuously varying j3- 
field. If we assume that each cell reaches local equilibrium very fast, i.e. the 
associated relaxation time is short, then in the long-term run the stationary 
probability distributions p{E) arise as the following mixture of Boltzmann 
factors: 

P{E) = f{(5)^^p{E)e-P^d(5 (2) 

Without restriction of generality, we may absorbe the factor 1/Z{I3) into the 
function i.e. define f{(5) — f{(5)/Z{l3) and rename f ^ f. Also, for 

reasons of simplicity we may just assume p{E) = 1, keeping in mind that the 
most general case may correspond to a different density of states. The result 
is an effective distribution 

POD 

p{E) ~ / /(/?)e-^^d/? (3) 

given essentially by the Laplace transform of 

The simplest dynamical example of a superstatistical system is a Brown- 
ian particle of mass m moving through a changing environment in d dimen- 
sions. For its velocity v one has the local Langevin equation 

ij ^ --fv + aL{t) (4) 
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{L{t): (i- dimensional Gaussian white noise, 7: friction constant, a: strength 
of noise) which becomes superstatistical if the parameter f3 := is re- 
garded as a random variable as well. In a fluctuating environment (3 may 
vary from cell to cell on a large spatio-temporal scale. Of course, for this ex- 
ample the energy E is just kinetic enery E = ^mv^. In each cell of constant 
P the local stationary velocity distribution is Gaussian, 

Pim=(j^y'e--^^-^, (5) 

provided the relaxation time 7"^ is small enough as compared to the changes 
of (3. The marginal distribution describing the long-time behaviour of the 
particle 

POO 

Piv) = / f{P)p{mdP (6) 



exhibits fat tails. The large-|w| tails of the distribution ^ depend on the be- 
haviour of f{P) for /? — i> [F. Different superstatistical models correspond- 
ing to different The function / is determined by the environmental 
dynamics of the nonequilibrium system under consideration. 

Consider, for example, a simple case where there are n independent Gaus- 
sian random variables Xi, . . . , X„ underlying the dynamics of (3 in an additive 
way. P needs to be positive and a positive /? is obtained by squaring these 
Gaussian random variables. The probability distribution of a random vari- 
able that is the sum of squared Gaussian random variables (3 = ^17=1 -^f 
is well known in statistics: It is the ^^-distribution of degree n, i.e. the 
probability density is given by 

where jSo is the average of (3. In this case the integration over the fluctuating 
f3 can be explicitly done and one obtains 



P{v) z (8) 

{1 + l3{q - l)lmiP)~ 



with 
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and 



Wo 



(10) 



2-{q-l)d' 



These types of distributions, so-called g-exponentials, are relevant for nonex- 
tensive statistical mechanics [H 121 El Ej . 

For other systems the random variable (3 may be generated by multi- 
plicative random processes. In such cases one typically ends up with a log- 
normally distributed (3, i.e. the probability density is given by 



where and s are parameters. The corresponding superstatistics is relevant 
for turbulent systems |2Hl 1211 jHOj ■ The integration over [3 cannot be done 
analytically in this case. 

3 Generalization 

Let us now generalize the additive case. For more general nonequilibrium 
situations, it may not (3 itself that is x^-distributed, but some suitable power 
of j3, say y := where is a parameter. What the relevant power r] 
is may depend on the physical nature of the nonequilibrium system under 
consideration and the physical meaning of 13. rj can depend on the spatial 
dimensions of the experiment under consideration, its boundary conditions, 
the dissipation mechanism, the nature of the driving forces, etc. 

To construct a microscopic model for y = j3^, we may assume that 
there are many (nearly) independent microscopic random variables j = 
1, . . . , J, contributing to y in an additive way. For large J their rescaled sum 

X]/=i ^ill approach a Gaussian random variable Xi due to the Cen- 
tral Limit Theorem. In total, there can be many different random variables 
consisting of microscopic random variables, i.e., we have n Gaussian random 
variables Xi, . . . ,Xn due to various relevant degrees of freedom in the sys- 
tem, y must be positive. A positive value is obtained by squaring the X^. 
By construction, the resulting y = ^"^^^ Xf is x^-distributed with degree n, 
i.e. one has the probability density 




(11) 




(12) 
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where yo is the average of y. 

A priori, all kinds of values of rj are possible. For example, rj = —1 
means that the temperature is the relevant quantity that can be repre- 
sented as a sum of several squared Gaussian random variables. Another case 
may be that the local standard deviation of the Gaussians is x^-distributed, 
corresponding to the case rj = — |, and so on. Each choice of t] implies 
a different probability distribution of (3. The resulting probability density 
fpiP) can simply be obtained by employing a transformation from the x^- 
distributed random variable y = (3^ to the original random variable (3 via 
fy{y)dy = fp{(3)df5. One obtains 

where y^ = {y) = (/?'') is the average of y. As we shall see in the next 
section, for 77 > this generates asymptotic power laws for p{E), whereas 
the case 77 < generates stretched exponentials. Both cases can be physically 
relevant. 



4 Asymptotic behaviour 

Let us study the marginal probability distributions p{E) as given by eq. Q 
for our class of inverse temperature distributions ()13p . In particular, we are 
interested in the behaviour of p{E) for large E, i.e. the shape of the tails. 

Two important cases have to be distinguished, rj > and rj < 0. As 
has been shown in |Hj, the large-energy tails of p{E) are determined by the 
behaviour of the function for /? — > 0. For 77 > eq. (fT^ yields the 
small-/? asymptotics 

MP) - P'^-\ (14) 

since the exponential term exp{—^i3^) in eq. (fT^ just approaches 1. Ac- 
cording to the general results of jH] this implies 

p{E) ~ E-^t (15) 

for E ^ 00. This means one obtains an asymptotic power-law decay in 
E with the exponent t]^. For example, if the energy is just kinetic energy, 
E = |mf^, then p{v) decays with tails proportional to |f |~''". This reminds 
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us of Tsallis statistics and nonextensive statistical mechanics PQ 121 El El, with 
an entropic index q given by 



2 2 

= rfn <^ q = l-\ . (16) 

g — 1 rjn 

However, only 77 = 1 yields exact Tsallis statistics, whereas the p{E) obtained 
for other 77 > behave slightly different for finite values of E. Asymptotically, 
however, only the product of the two parameters rj and n is relevant. The 
same asymptotic power law can be achieved by e.g. doubling the number 
of degrees of freedom n if at the same time the exponent rj is reduced by a 
factor 1/2. 

The case 77 < is very different. Here for small (3 the behaviour of the 
function is dominated by the exponential term 

fm-e-'^P\ (17) 

According to the general results of ^ , this implies the large energy behaviour 

p(E)~e«^^, (18) 
where a is a negative constant^ depending on c = 2^ and 77: 

1 + ^ 

(19) 



(cIt^I)"- 



We see that p{E) now asymptotically decays as a stretched expontial. It is 
interesting to see that the degrees of freedom n do not influence the exponent 
of the stretched exponential, they only influence the proportionality constant 
a. 



5 Applications 

For granular gases, stretched exponential tails (and tales!) are known to 
play an important role ^1 031. For example, Rouyer and Menon 
performed an experiment where they measured the probability distribution 
of velocities in a granular system confined to a vertical plane and driven 

"'^ There is a misprint in eq. (29) of 'E, which is hereby corrected. 
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by strong vertical vibration. The system reaches a stationary state with a 
probabiUty density of horizontal velocities given by 

~e-^l"l", (20) 

where a is measured to be a = 1.55 ± 0.1. This result is robust for all 
frequencies and amplitudes. The theoretical explanation of this behaviour is 
still unclear. 

We may assume that a superstatistics of the type described in this paper 
is relevant. That means we assume that some suitable power y = (3"^ of the 
inverse granular temperature P is x^-distributed due to large-scale spatio- 
temporal inhomogenities of granular temperature. We obtain agreement with 
the experimentally observed stretched exponential tails if i] is given by 

2 — '— = a^ri = ^-3.4. (21) 

Of course a complete theory would have to explain why in this experiment 
granular temperature raised to a power of approximately 3.4 is the relevant 
X^-distributed observable. Currently this is out of reach. The value rj = —3.4 
seems not to be universal. Molecular dynamics simulations [IHl of dilute 
granular gases driven by vertically oscillating plates yield exponents a in the 
range 1.2 < a < 1.6, which is equivalent to —4.0 < rj < —2.3. 

Stretched exponentials are sometimes also used to fit data in hydrody- 
namic turbulence, despite a lack of theory for this. For example, Bodenschatz 
et al. |2S| present a fit of the measured probability distributions p{a) of ac- 
celeration a of a Lagrangian tracer particle which asymptotically decays as 
a stretched exponential, 

pia) ~ e-'=l'^l°'. (22) 

This behaviour, interpreted in terms of a superstatistics, would correspond to 
1] = —1/4. This would mean that the square root of the standard deviation 
of the local Gaussians would be the relevant x^-distributed random variable. 

Overall, one has to be careful with over-ambitious fits of experimental 
data using stretched exponentials. Power laws as given by g-exponentials, 
lognormal superstatistics and stretched exponentials can all yield very sim- 
ilar looking fits of fat-tailed distributions if the largest accessible energy E 
lies within 5-10 standard deviations. In typical experimental situations the 
asymptotics of the tails is often not reached. So there can be many competing 
models to explain the data. 
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